The Fermi-polaron in two dimensions: Importance of the two-body bound state 
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We investigate a single impurity interacting with a free two-dimensional atomic Fermi gas. The 
interaction between the impurity and the gas is characterized by an arbitrary attractive short- 
range potential, which, in two dimensions, always admits a two-particle bound state. We provide 
analytical expressions for the energy and the effective mass of the dressed impurity by including 
the two-body bound state, which is crucial for strong interactions, in the integral equation for the 
effective interaction. Using the same method we give also the results for the polaron parameters in 
one and three dimensions finding a pretty good agreement with previous known results. Thus our 
relations can be used as a simple way to estimate the polaron parameters once the two-body bound 
state of the interaction potential is known. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.-d 
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I. INTRODUCTION 

In the last years, thanks to the experimental capabil- 
ity of tuning the relative population and the interaction 
strength in atomic gases of different species, highly im- 
balanced gases have been extensively studied. In par- 
ticular a lot of theoretical and experimental work has 
been devoted to atomic Fermi gases in two different hy- 
perfine states (see, e.g., the recent review [l| and refer- 
ences therein). The building block is the solution of the 
limiting case of a single impurity atom interacting via 
a short range potential with an ideal atomic Fermi gas. 
Such a problem is not only relevant in the field of ultra- 
cold gases, since it is related to the more general one, 
the so-called impurity problem, which is present also in 
other area of physics. In our case the dressed impurity is 
called Fermi-polaron (or polaron) in analogy with elec- 
trons dressed by the bosonic (phonon) bath in a crystal. 

Important quantities characterizing the polaron are: 
(i) its chemical potential, also called interaction energy 
or binding energy, i.e., the (negative) energy difference of 
the ground state with and without the impurity atom at 
rest; (h) its effective mass, i.e., the dressed parabohc dis- 
persion relation at low impurity's momentum. In three 
dimensions (3D) these parameters have been calculated 
in different ways, e.g., by means of variational ansatz 
d, 0], Monte-Carlo methods functional renormal- 

isation group @ and they have been measured experi- 
mentally 0- 

The variational approach is known to give reasonable 
results also in the one-dimensional (ID) case 0, since 
it can be compared with the exact solution found by 
McGuire 0. Very recently the same approach has been 
used to study the two dimensional (2D) case, where it 
has been shown that in 2D its use can be questionable at 
least if using the same approximations as in 3D 10, 11]. 

In two-dimensions, as well as in one dimension, an at- 
tractive interaction always allows for a two-body bound 
state. In the present work we solve the impurity prob- 
lem including such a bound state explicitly in the integral 
equation for the effective interaction of the impurity with 



the Fermi gas. Within a number of approximations we 
can provide analytical expressions for the polaron param- 
eters which agree quite well with the known results in one 
and three dimensions. In 2D we find an expression for the 
energy which interpolates between the correct/expected 
limiting values in the weakly and the strongly interacting 
regime. Thanks to the recent experimental advances in 
realizing one- and two-dimensional strongly interacting 
Fermi gases [l3| , the impurity problem in reduced di- 
mensionality has become relevant in the context of ultra- 
cold gases. 

In the next section we introduce the formalism and 
give the result for the energy and the effective mass of 
the two-dimensional Fermi-polaron problem. In Sec. IIIII 
we apply the method to the one- and three-dimensional 
case. 



II. FERMI-POLARON IN 2D 

It is known that Brueckner-Hartree-Fock theory, when 
applied to the Fermi-polaron problem in three dimension 
gives reasonable results [iJl (see also Sec. IIIIBp . The 
basic equation from this theory is the Bethe-Goldstone 
integral equation for the reaction matrix [l^ , also called 
effective interaction, (e.g. for the 2D electron gas flGf). 
The Bethe-Goldstone integral equation for the effective 
interaction between a particle in the bath with momen- 
tum ki and the impurity atom with momentum k2 can 
be written as 
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In Eq. ([T]) q is a transfer momentum, ^(q) is the 
Fourier transform of the two-particle interaction poten- 
tial, Tik the Fermi distribution function at zero temper- 
ature. The interaction energy or correlation energy fol- 
lows then from the mean value of the effective interaction 
ep = (5(ki,k2,q)) 
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For k2 = one gets the rest correlation energy 6° of 
the polaron, and by expanding about this solution in fc| 
one get its effective mass m* as usual by the relation 
E = e° + fc|/2m*. We remind that in Eq. ([T]) for the 
effective interaction only ladder diagrams are summed 
and the Fermi sea limits the momenta in the intermediate 
states. 

We consider an interaction characterized by an at- 
tractive short-range potential of arbitrary shape. In 
three dimensions this potential can be approximated by 
a delta function, and Eq. ([1]) coincides with the self- 
consistent equation obtained via single particle-hole vari- 
ational ansatz 0, In two dimensions it is not clear, 
whether one can use a delta function pseudo-potential, 
(see e.g. hence, a solution of Eq. ([T|) obtained in 

the same way as in 3D is questionable. As usual in order 
to treat properly the 2-body problem we write Eq. ([!]) 
by expressin g V in terms of the two-particle scattering 
amplitude / [1^ 



S((ki,k2,q) = /(ki - k2,q) - 
dk /(ki - k2,q- k) (-riki+k) 
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where /(ki — k2,q) is the off-shell scattering amplitude 
[l9j . Note, that this equation is already renormalized 
with repect to ultraviolet divergencies. 

For short-range potentials the exchange momentum 
q is small and the main contribution from the integral 
comes from small values of k, thus we can approximate 
the off-shell scattering amplitude by the on-shell scatter- 
ing amplitude /(ki — k2). Then the effective interaction 
does not depend on the exchange momentum and for the 
impurity at rest (fc2 = 0) Eq. ^ reduces to 



5(ki) 



/(ki) 



dp 



p<ki 



(2^)2p(p-ki) 
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We assume that the finite range R of the attractive 
interaction potential is the shortest length scale in the 
system. In particular for kpR ^ 1 the s-wave scattering 
amplitude reads (see e.g. 



f{k)-^ = -[ln(fcV™|eb|) - i7r]m/47r 



(4) 



where e^, is the binding energy of the two-particle bound 
state, which in 2D (differently from 3D) is always present 

Solving Eq.® for the weakly interacting case |eb| <C 
ep, where ep = h'^kp/(2m) is the Fermi energy, we find 
the interaction g{ki) = — 47r/[mln(2e^/|eti|)] from which 
we get the mean field energy 6° = — 2eF/[?Tiln(2e_F/|ef,|)], 
which is obviously in agreement with the result found in 
the weakly interacting regime using the single particle- 
hole variational ansatz [10[. 

Until now, deriving Eq. ^ from the Bethe-Goldstone 
Equation, we have neglected the two-particle bound 



state. In order to take it into account we go back to Eq. 
([2]). We rewrite the initial energy of excitation processes 

k'^ k^ p2 

appearmg m the denommator as 2^ + 2m ~ "m~'~im' 
with relative momentum k^ and centor-of-mass momen- 
tum P. Further, we remind that the scattering ampli- 
tude /(kr) depends on the relative momentum only. It 
is well known that in the case where the majority par- 
ticle with ki and the impurity with k2 form a two-body 
bound state, their relative momentum kr is purely imag- 
inary with — = £{,< being the binding energy. For 
the calculation of the interaction energy we assume the 
impurity to be at rest and P = 0. Then the effective 
interaction in ladder approximation obeys instead of Eq. 

m 
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Let us notice that a very similar equation is found in [23| 
for the vertex function in the presence of a two parti- 
cle bound-state, where the molecular propagator is ex- 
pressed by the two-particle scattering amplitude at the 
vacuum energy of the molecule. 

At the momentum corresponding to the bound state 

the scattering amplitude has a pole / {kr — ^/rncb) = 
and the interaction energy is given by 



ng{eb) 



and it is shown in Fig. [T] 
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Polaron energy and effective mass in 2D 
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FIG. 1: Main: polaron energy as a function of the two-body 
binding energy et as given by Eq. ^ (solid-red line). For 
completeness we report also the weakly interacting (dashed- 
green line) and the strongly interacting (dotted-black line) 
results (see text). Inset: ratio m*/m between the effective 
and the bare mass as given by Eq. 



In the limit of weak interactions, i.e., |eb| <C ep, the po- 
laron energy reduces to « — 26^/ ln(2ei?/|eb|), which 
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coincides with the one obtained from Eq. ([3]) or using 
the variational ansatz [loj , where the two-particle bound 
state is not included. In the opposite limit, |e;,| cp, 
Eq. 1^ yields ep = —\eb\ — ep + o(e^-/|efc|). This is the 
expected result, because in this regime one atom of the 
majority is strongly bound to the impurity with — |eb| 
and thus it has to be removed from the Fermi-sea lead- 
ing to the first correction —ep- Thus Eq. ^ smoothly 
interpolates between these two limits and it provides a 
good approximation for the interaction energy of the 2D 
polaron in all regimes. Our results are in good agree- 
ment with recent preliminary Monte-Carlo calculations 
[11]. We remind that Eq. © is valid for all attractive 
potentials with s-wave scattering amplitude of logarith- 
mic form Eq. ([4]). 



A. Effective mass 

In this section we study the effect of the interaction 
on the motion of the impurity. As already mentioned 
its energy can be expanded for small momentum k2 as 
E = ep + fc|/2TO*, where we define the effective mass of 
the impurity as 

1 

TO* 



1 

m 



1 &p(k2) 
mdm/2) 



\k2=0 



(7) 



and ep(fc2) — ng(eh,k2), with the effective interaction 
g calculated again including the two-body bound state. 
In particular the initial energy of excitation processes is 
— \^b \ + ^1/4™ and instead of Eq. (O we obtain 



1 



5(ef.,k2) 



1 



f{kr) Jo (27r)2 |et,| + e(fc2,fc) 



(8) 



wher 



k2 k cos 



-f- — . The ratio between 



ere e(fc2,fc) = 4^ + 
the bare and the effective mass of the impurity atom in 
two dimensions within our approximation reads 



m 
m* 




(9) 



As shown in Fig. [l]the effective mass to* obtained from 
the previous equation has the expected behavior: for 
small interactions it is close to the bare mass value to 
and for large interactions it approaches the molecular 
mass value 2to. 



III. FERMI-POLARON PARAMETERS IN ID 
AND 3D 

For a two-dimensional system our approach seems to 
give quite reasonable results and provides analytical ex- 
pressions for the polaron's parameters. In the present 
section we apply our approach to the one- and three- 
dimensional case. The simple expressions we find are in 
reasonable agreement with the known results. 



A. The Fermi-polaron in one dimension 

In one dimension the Fermi-polaron problem admits 
an exact solution and the interaction energy reads 



p 

2eF 



y + ^y 



(1 + y^) arctan(y) 



(10) 



where y = \/\eb\/ (2e_F). Again e;, is the binding energy 
of the lowest two-body bound state. 

When applied to ID Eq. ^ gives for the polaron 
energy 



p 

2eF 



y 



arctan(i) 



(11) 



which we compare against the exact result Eq. ([10)) in 
Fig. [2j The agreement looks pretty good, although in the 

Polaron energy and effective mass in ID 




FIG. 2: Interaction energy (main panel) and effective mass 
(inset) as a function of the two-body binding energy et in ID. 
The approximate result Eq. (solid-red line) is cornpared 
with the exact expression Eq. (|10|) given by McGuire 



strongly interacting case, we get — |e;,| — 2/3eF instead of 
— |eb| — ep- However our results are closer to the exact 
solution than the one obtained in with the single particle- 
hole variational ansatz (see e.g., Q). 



B. The Fermi-polaron in three dimensions. 

In a three-dimensional geometry when no two-body 
bound state is present, i.e., for negative s-wave scat- 
tering length (a < 0), we can use Eq. ^ - the usual 
Brueckner-Hartree-Fock theory - with scattering ampli- 
tude 47r/(ki)~^/TO = a""'^ -|-i|ki|. The main contribution 
to the effective interaction is g{ki = 0) and thus we can 
write an approximated expression for the polaron energy 
as 



2 kpa 



(12) 
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Polaron energy in 3D 
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FIG. 3: Interaction energy as a function of the inverse 3D 
scattering length l/(feFffl) (red line) in comparison with the 
results obtained from Monte-Carlo calculations 0] in the 
polaron regime and at unitarity (black diamonds). In the 
molecule regime (|ei,| 3> tp) we compare the expected result 
ei„t = —\et\ — (.F (blue line). Inset: zoom on the negative 
axis. 



Note, that adding self-consistency to the Bethe- 
Goldstone equation by changing the initial energy of exci- 



2m 2m ^ 



increases 



tation processes 2m~^ Im ^ 
the accuracy of the results in the polaron regime and 
yields the same equations as in 0]. 

In the molecular regime the potential admits a two- 
body bound-state (a > 0) with binding energy e;, = 
-l/(ma2). From the 3D version of Eq. ([S]) we obtain 
the interaction energy of the impurity, which reads 



p 



(13) 



3 1- 



■ arctan I 



with Y^lebl/ (2ei?) — l/{kpa). For large binding energy 
Icbl > ef (or kpa < 1), one gets e°„t = -|eh| - 66^/5, 
which is larger than the expected result by —1/5 e^- In 
Fig. [3] we show the comparison between the previous 
simple expression and the results obtained from Monte- 
Carlo calculations [H . Although the agreement is good it 
is worse, as already mentioned, than the results obtained 
from the variational approach in both the polaron [2] and 
the molecular regime [3| as well as the results obtained 
from the functional renormalization group [S] . 



where y = ^y\eb\/{2ep) is defined as in Eq. ([10]). Again 
the agreement is reasonable and better than the single 
particle-hole variational ansatz 0, . 

In 3D the situation is more involved since a maximum 
(quite larger than 2™) for the effective mass has been 
found when the nature of the impurity changes from a 
fermionic quasi-particle to a bosonic quasi-particle 
In our approximation this maximum cannot be found. In 
order to find this maximum one has to take three-particle 
scatteirng into account which is beyond the scope of this 
paper. 



IV. CONCLUSION 



In conclusion, we have investigated the problem of an 
impurity atom interacting with a non-interacting Fermi- 
gas in a two-dimensional geometry. We consider a short 
range, attractive potential, which implies the presence 
of a two-bound state for any interaction strength. We 
have calculated the interaction energy and the effective 
mass of the impurity by including the bound-state in the 
Bethe-Goldstone integral equation for the effective inter- 
action. We were able to obtain simple analytical expres- 
sions which give reasonable results in 2D, interpolating 
between the weakly and the strongly interacting regime 
(see Fig. [Ij. Moreover when applied to the three- and 
one-dimensional case our polaron parameters well com- 
pare with most of the known results (see Sec. IIII|) . Thus 
our expressions can be used to estimate the dressed im- 
purity's parameters in a simple way once the two-body 
bound state is known. Our analysis shows how important 
the two-body bound state is for the polaron problem. 

Finally let us stress that we do not address the question 
of whether there exists a polaron-to-molecule transition 
in two-dimension as debated in [l^, [HI ■ The problem is 
still open and it could happen that the system behaves 
similarly to the one-dimensional case where there is not 
such a transition. Such a question is clearly relevant for 
the possible low-temperature phases of a two-dimensional 
polarized Fermi gas, whose balanced version has been 
recently experimentally realized [1^. 



Remarks on the effective mass in ID and 3D 
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The effective mass of the impurity is given by Eq. ^ 
in any dimension. 

In the inset of Fig. [2] we compare our result (Eq. [9]) 
with the exact one found by McGuire in [§| which reads 

m* _ (l + I arctan y)^ 
1 + 1 (arctan y-h^ 
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